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Abstract. Some important problems (e.g., in optimal transport and optimal control) have 
a relaxed (or weak) formulation in a space of appropriate measures which is much easier 
to solve. However, an optimal solution y of the latter solves the former if and only if 
the measure y is supported on a “trajectory” {(t,x(t)): t £ [0, T]} for some measurable 
function x[t). We provide necessary and sufficient conditions on moments (7y) of a measure 
dy(x, t) on [0, l] 2 to ensure that y is supported on a trajectory {(t, x(t)): t £ [0,1]}. Those 
conditions are stated in terms of Legendre-Fourier coefficients f/ = (fj(i)) associated with 
some functions fj: [0,1] —1 K, j = 1,..., where each f/ is obtained from the moments 7^, 
i = 0,1,..., of y. 
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1 Introduction 

This paper is in the line of research concerned with the following issue: which type and how 
much of information on the support of a measure can be extracted from its moments (a research 
issue outlined in a Problem session at the 2013 Oberwolfach meeting on Structured Function 
Systems and Applications [2]). In particular, a highly desirable result is to obtain necessary 
and/or sufficient conditions on moments of a given measure to ensure that its support has 
certain geometric properties. For instance there is a vast literature on the old and classical 
L-moment problem, which asks for moment conditions to ensure that the underlying measure y 
is absolutely continuous with respect to some reference measure 17 and with a density in L 00 (i/). 
See, for instance, [3, 10, 11], more recently [7], and the many references therein. 

Here we are interested in a problem that is somehow “orthogonal” to the L-moment problem. 
Namely, we consider the following generic problem: Let dy(x, t ) be a probability measure on 
[0,1] x [0,1]. Provide necessary and/or sufficient conditions on the moments of y to ensure 
that y is singular with respect to the Lebesgue measure d(x,t ) on [0,1] 2 - In fact, and more 
precisely, suppose that: 

• one knows all moments 7j(j) = f x l t ] dy(x, t ), i,j = 0,1 ,..., of the measure y, and 

• the marginal of y with respect to the “t” variable is the Lebesgue measure dt on [0,1]. 

Then provide necessary and/or sufficient conditions on the moments (7 i(j)) of y to ensure 
that y is supported on a trajectory {(t,x(t)): t £ [0,1]} C [0,1] 2 , for some measurable function 
x : [0,1] -A [0,1]. 

*This paper is a contribution to the Special Issue on Orthogonal Polynomials, Special Functions and Applica¬ 
tions. The full collection is available at http://www.emis.de/journals/SIGMA/OPSFA2015.html 
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In contrast to the L-moment problem, and to the best of our knowledge, the above problem 
stated in this form has not received a lot of attention in the past even though it is crucial in 
some important applications (two of them having motivated our interest). 

Motivation. In addition of being of independent interest, this investigation is motivated by 
at least two important applications: 

The mass transfer (or optimal transport ) problem. In the weak (or relaxed) Monge- 
Kantorovich formulation of the mass transport problem originally stated by Monge, one searches 
for a measure dp(x, t ) with prescribed marginals v x and Uf, and which minimizes some cost func¬ 
tional f c(x, t ) dn(x , t). However in the original Monge formulation, ultimately one would like to 
obtain an optimal solution fi* of the form dfi*(x,t ) = 5 X ^ dvfit) for some measurable function 
1 1 -)- x(t) (the transportation plan) and a crucial issue is to provide conditions for this to happen 1 . 
For more details the interested reader is referred, e.g., to [13, pp. 1-5] and [9]. There exist some 
characterizations of the support of an optimal measure for the weak formulation. For instance, 
c-cyclical monotonicity relates optimality with the support of solutions, and more recently [ 1 ] 
have shown in the (more general) context of the generalized moment problem that under some 
weak conditions the support of optimal solutions is finitely minimal / c-monotone. (As defined 
in [ 1 ] a set T is called finitely minimal / c-monotone if each finite measure a concentrated on 
finitely many atoms of T is cost minimizing among its competitors; in the optimal transport 
context, a competitor of a is any finite measure a' with same marginals as a.) For more details 
the interested reader is referred to [1] and the references therein. But such a characterization 
does not say when this support is a trajectory. 

- Deterministic optimal control. Using the concept of occupation measures, a weak formu¬ 
lation of deterministic optimal control problems replaces the original control problem with an 
infinite-dimensional optimization problem V on a space of appropriate (occupation) measures 
on a Borel space x x [0,1] with Sfi C M n , % C M m . For more details the interested reader 
is referred, e.g., to [8, 14], and the many references therein. An important issue is to provide 
conditions on the problem data under which the optimal value of the relaxed problem V is the 
same as that of the original problem; see, e.g., [14]. Again this is the case if some optimal solu¬ 
tion (i* (or every element of a minimizing sequence) of the relaxed problem is such that every 
marginal y* of fi* with respect to ( Xj,t ), j = 1,..., n, and every marginal pf of fi* with respect 
to ( U£,t ), £ = 1 ,m, is supported on a trajectory {(t,Xj(t)): t G [0,1]} and on a trajectory 
{(t, ufit )): t G [0,1]} for some measurable functions t H > Xj(t) and t (->• ufit) on [0,1]. 

Contribution. Of course there is a particular case where one may conclude that p is singular 
with respect to the Lebesgue measure on [0,1] 2 . If there is a polynomial p G M[x,f] of degree 
say d, such that its vector of coefficients p is in the kernel of the moment matrix M s (where 
M s [(*,j), (k,£)\ = 7 i + k,j+e, i + j, k + £ < s, with d < s), then p is supported on the variety 
{(x,f) G [0,1] 2 : p(x,t ) = 0} and therefore is singular with respect to the Lebesgue measure 
on [0, l] 2 . But it may happen that p(x, t) = p(y, t) = 0 for some t and some x fi y and so even in 
this case additional conditions are needed to ensure existence of a trajectory {(£, x(t)): t G [0,1]}. 

We provide a set of explicit necessary and sufficient conditions on the moments qq = ( 7 i(j)) 
which state that for every fixed i, the moments 7 i(j), j = 0 , 1 ,..., are limits of certain i-powers 
of the moments 7 1 . 

More precisely, an explicit linear transformation of the infinite vector / ~/ 1 is the vector 
of (shifted) Legendre-Fourier coefficients associated with the function 1 1 - 7 - x(t). Then the con¬ 
ditions state that for each fixed i = 2,3,..., the vector A 7 j should be the vector of (shifted) 
Legendre-Fourier coefficients associated with the function 1 1 —> x(t) 1 , which in turn are express¬ 
ible in terms of limits of “i-powers” of coefficients of A^ 1 . 

At last but not least, it should be noted that all results of this paper are easily transposed to 
the multi-dimensional case of a measure dp(x, t) on [0, l] n X [0, 1] and supported on a trajectory 

1 Hcre 5, denote the Dirac measure at the point z. 
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{(f,x(t)): t € [0,1]} C [0,1]" +1 for some measurable mapping x: [0,1] —> [0,1]". Indeed by 
proceeding coordinate-wise for each function t i->- Xi(t), i = 1, ..., n, one is reduced to the case 
[0, l] 2 investigated here. 


2 Notation, definitions and preliminary results 

2.1 Notation and definitions 

Given a Borel probability measure p on [0,1] 2 , define 

i,j = 0,1,..., 


7 i{j)= / t J x l dp(x,t), 

J[ o,i ] 2 


( 2 . 1 ) 


and for every fixed i E N, denote by 7 , : the vector of moments (7 i(j)), j = 0,1,.. ■■ 

Let (Jfj), j = 0,1,..., be the family of orthonormal polynomials with respect to the Lebesgue 
measure on [0,1]. They can be deduced from the Legendre polynomials 2 via the change of 
variable t' = (2t—1); the (Jzf)) are called the shifted Legendre polynomials. The polynomials (Jzf)) 
can be also computed exactly from the moments 7 0 of the Lebesgue measure dt on [0,1] by 
computing some determinants of modified Hankel moment matrices. For instance, Jz?o = 1, and 

Jz?i(t) = odet = aft — 1/2) with a 2 (l/3 — 1/2 + 1/4) = 1, 

i.e., a = \/l 2 and (t) = 2\/3 1 — y/3, and 
( 1 1/2 l/3\ 

2 z? 2 (t) = 6 det 1/2 1/3 1/4 = b[t 2 /12 - t/12 + 1/72], 

V 1 t #) 

with b > 0 such that /q 1 Jz? 2 (t) 2 dt = 1. See, e.g., [4] and [6]. 

The Lebesgue space L 2 QO, 1]) := {/: [0,1] —> M, f 2 dx < 00 }, equipped with the scalar 
product (/, g) = fgdx and the associated norm || ■ ||, is a Hilbert space and the polynomials 
are dense in L 2 ([ 0 ,1]). In particular the family (Jz fj), j = 0,1,..., form an orthonormal basis 
of L 2 ([0,1]). Let i 2 denotes the space of square-summable sequences with norm also denoted 

by II • II- 

Finally, let ||/||oo := ess sup \f(x)\, and similarly, for p E M[x] let ||p||oo := sup |p(x)|. 

is [ 0 , 1 ] ®e[o,i] 

Then ^([0,1]) := {/: ||/||oo < 00 }. 


2.2 Some preliminary results 

We next state some useful auxiliary results, some of them being standard in Real Analysis. 
Proposition 2.1. Let 1 1 —> f(t ) be an element o/L 2 ([0, 1]) and define f = (f(j)) by 

f(j) := f &j(t)f(t)dt, j = 0,1,.... 

■Jo 

Then one has 


( := lim J = / in L 2 ([0,1]), 


n —>00 

3=0 \ j =0 


( 2 . 2 ) 


and this decomposition is unique. Moreover f G f 2 and 




2 The Legendre polynomials are orthonormal w.r.t. to the Lebesgue measure on [—1,1]. 
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When the interval is [—1,1] (instead of [0,1] here) (2.2) is called the Legendre (or Legendre- 
Fourier) series expansion of the function / and f = (f (/)) is called the vector of Legendre-Fourier 
coefficients. 

The notation f k stands for the function t (-»• f(t) k , k £ N. If f k £ L 2 GO,1]) we denote by 
ffc = (ffc(j)) £ its (shifted) Legendre-Fourier coefficients so that ||/ fc || = |f^ || (where again 
the latter norm is that of l 2 ). Notice that we also have: 

Proposition 2.2. Let f,f k ,g G L 2 GO,1]) with k G N. Then g = f k if and only if g = f k . 

This follows from the uniqueness of the decomposition in the basis (Jz fj). 

We also have the following helpful results: 


Lemma 2.3. Let f G L 2 GO, 1]) with Legendre-Fourier coefficients f = (f(j))> and let ( p n ) C M[t] 
he a sequence of polynomials such that \\p n — f\\ —> 0 as n —> 00 . 

If Pn = (PnU )) denotes the (shifted ) Legendre-Fourier coefficients of p n , for all n = 1, 2 ,..., 
then || p n — f || — > 0 in l 2 as n 00. 

Proof. As II p n - /II 2 = fo(p n ~ f) 2 dt and with d n = deg (p n ), 




d n /»1 /* 1 

^(Pnti)) 2 / dt+2^p n {j)p n {k) / SfjSfkdt 

i = 0 -- fc<J ^-..-- 


=1 

fi 


=0 


2 ^Pn{j) / &jfdt + 


j =0 


m 


=l|f||=E 7 f(i ) 2 


> Y^Pntif ~ 2 PnU)IU) + f 0') 2 = £(Pn(j) “ f{j)Y 


3=0 


3=0 


and the result follows because ||p n — /|| —> 0 as n —> 00. ■ 

Lemma 2 . 4 . Let f G LqoQO, 1 ]) (hence f k G L2QO, 1 ]) for every k G N). If a sequence (p n ) C 
M[x] is such that sup n ||p n ||oo < 00 and || p n — /|| —> 0 as n 00 , then for every k G N, 



In addition, if p k n denotes the ( shifted ) Legendre-Fourier coefficients of p k , n = 1,..., then 
|| p^ — ffc|| —>■ 0 in i 2 as n ^ 00. 


Proof. Let M > max 


3 ,sup n ||p n ||oo] and fix k G N, 


Pn~f k \\ 2 = I {Pn f k ) 2 dx = / {p n -f) 2 [^2p k -^f ] dx 


pk\ 2 


fk -1 


£=0 


< / (Pn -/) 2 dx<{kM k - 1 ) 2 / ( Pn -f) 2 dx 


\£=1 


= (kM k 1 ) 2 ||p n -/|| 2 (-t 0 as n 4 00 .) 


Then the last statement follows from Lemma 2.3. 
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Definition 2.5. Let / G -^([0,1]) with Legendre-Fourier coefficients f. For every k,n G N, 
define the polynomial fn' > G M[x] and the vector G M fcn+1 by 

nk 

t tf’w := I £famm I = £>W)-s5(0- 

D =0 / j =o 

Observe that each entry fn^O), j = 0,..., n/c, is a degree-/c form of the first n + 1 Legendre- 
Fourier coefficients of f. Completing with zeros, consider f //'"" 1 to be an element of l 2 and if f 
converges in i 2 as n —>■ oo, call G i 2 its limit. 

The limit f 1 ^ can also be denoted f * • • • * f , the limit of the k times “^-product” in i 2 of the 
vector f G f 2 by itself. Equivalently one may write f® = since fn\t ) = fii~ l \t)fn\t) 


for all t G [0,1], and fiO ^ —» f fc 1 as n —> oo, as well as f, 


(i) 


f. 


Lemma 2.6. Let f G Loo([0, 1]), hence f k G ^([0, 1) with ( shifted ) Legendre-Fourier coeffi¬ 
cients ffc G £ 2 for every k G N, and assume that 


SU P 11 Jn 11 oo 
n 


= sup||X»^ 


3 l|oo 


< OO. 


3=0 


Then = f * • • • * f (k times ) /or every k = 1,2,..., meaning that for every fixed k G N, 


lim 

n —/oo 


£f(im] 

0=0 


= lim 

n—foo 


= lim 

n —/oo 


Ain 




1=0 


I> 0 ‘)- 2 W* 

3=0 


= 0 . 


Equivalently, f& = f/-_i * f for every k = 2, 3,.... 

Proof. The result is a direct consequence of Lemmas 2.3 and 2.4 with = fn ' 1 (and the 
definition of the limit product” in Definition 2.5). 


3 Main result 

Assume that we are given all moments of a nonnegative measure dn(x, t ) on a box [a, b] x [c, d] C 
M 2 . After a re-scaling of its moments we may and will assume that ^ is a probability measure 
supported on [ 0 , l ] 2 with associated moments 

7 i(j)= xH J dfi(x,t), i,j = 0 , 1 ,.... 

J [o,i] 2 

We further assume that the marginal measure m with respect to the variable t, is the Lebesgue 
measure on [ 0 , 1 ], that is, 70 (j ) = l/(j + 1 ), j = 0 , 1 ,.... 

A standard disintegration of the measure n yields 


7 *0') = / t j x l dfi(x,t) = / / x>(dx|t) eft, i,j = 0, 1 ,..., 

J [ 0 , 1] 2 40 \ J [ 0 , 1 ] / 


(3.1) 
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where the stochastic kernel ip(-\t) is the conditional probability on [0,1] given t E [0,1]. Observe 
that the measurable function fi in (3.1) is nonnegative and uniformly bounded by 1 because 
\x l \ < 1 on [0,1] for every i, and so f % E L oo ([0,1]) for every i = 1,.... 

The vector j i = ( 7 i{j)), j = 0,1,..is the vector of moments of the measure = fi(t)dt 

on [0,1], for every i = 1, 2,.... The (shifted) Legendre-Fourier vector of coefficients q L of fi are 
obtained easily from the (infinite) vector 7 via a triangular linear system. Indeed write 

j 

-2 'jit) = J2 A, k t k , Vf E [0,1], j = 0,1,..., 

k =0 


where A ?J > 0, or in compact matrix form 


'-2b' 


'1' 

-2i 

= A 

t 



t n 





(3.2) 


for some infinite lower triangular matrix A with all diagonal elements being strictly positive. 
Therefore 


A 7* = A 



T 


'-2b' 

t 

r 1 

-2i 


fi(t)dt = 


t n 

Jo 






fi(t) dt = q % . 


Suppose that the measure ji is supported on a trajectory {(t,x(t)): t E [0,1]} C [0, l ] 2 for 
some measurable (density) function x : [0,1] —> [0,1]. The measurable function t 1 —> x(t) is an 
element of Loo([0,1]) because H^Hoo < 1. Then by Proposition 2.1, 

OO 

x = in L 2 ([0,1]), (3.3) 

j =0 

where x = (x(j)) E l 2 is its vector of (shifted) Legendre-Fourier coefficients (with ||cc|| = ||x||). 
Similarly, for every k = 2,3,..., the function 1 1-7 x{t) k is in Loo([0,1]) and 

OO 

x k = ^2x k (j)^j in T 2 ([0,1]), 
i=o 


with vector of (shifted) Legendre-Fourier coefficients x k E £ 2 such that ||x^j| = ||*fc||. 

We also recall the notation xff' E M fcn+1 for the vector of coefficients in the basis (Jifj) of 
( n V 

the polynomial x(j)«S^ (t) , and when considered as an element of f 2 (by completing 

v=° ' J 

with zeros) denote by x^ E l 2 its limit when it exists. 

Theorem 3.1. Let ^ be a Borel probability measure on [0, l] 2 and let 7j(j), i,j = 0,1, ..., be 
the moments of (jl in (2.1). 
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(a) If jjL is supported on a trajectory {(t,x(t)): t G 


function t *->• xft) on [0,1] and if sup 


3=0 


0,1]} for some nonnegative measurable 
< oo, then 

OO 


A ^ /•. A. /V A 

Xi = X y J = X * • • • ★ X j — Xi -1 * X , 

^ V 

i times 


Vi = 2 , 3 ,_ 


(3.4) 


Equivalently 

= (A7 1 ) W = (A7 i _ 1 * A7 1 ) W , Vi = 2,3,..., (3.5) 

where A is f/ie non singular triangular matrix defined in (3.2). 

(6) Conversely, if (3.5) ZioWs then p is supported on a trajectory {(t, x(t)): t G [0,1]} for some 
measurable function 1 1 —> xft) on [0,1], and (3.4) also holds. 


Proof. The (a) part. As p is supported on [0, l] 2 one has Hxlloo < 1 and so the function 
t xft) 1 is in L 2 ([0,1]) for every i = 1, 2,.... So let t H > xft) be written as in (3.3). Consider 
the function 1 i-)- xft) 1 , for every fixed i E N, so that 


OO 


j=o 


in L 2 ([0,1]), 


where the (shifted) Legendre-Fourier vector of coefficients x % is obtained by Xi = A 1 7 ,-. But 
by Lemma 2.6, we also have 


x l = 


2*0')-^- 
J=o 


lim 

77,—yoo l f 


in L 2 ([0,1]), 


J=o 


with ||x*|| = ||xj|| and x t = A^q. In other words, = Xj or equivalently, x 1 ' 1 ' 1 = A.'y i = 
(A 7 1 )^, which is (3.4). 

We next prove the (b) part. By the disintegration (3.1) of the measure p, 


7*0') = f t J fi(t)dt, i,j = 0,1,... 

Jo 

for some nonnegative measurable functions /* G Loo([0,1]), i = 1, 2,.... 

As (3.5) holds one may conclude that q i = q[ 1 ' where q { is the (shifted) Legendre-Fourier 
vector of coefficients associated with /* G L 2 ([0,1]), i = 1,.... Hence by Proposition 2.2, 
fi(t) = fifty a.e. on [0,1], for every i = 1,2,.... That is, for every i = 1,2,..., there exists 
a Borel set Bi C [0,1] with Lebesgue measure zero such that fift) = fift) 1 for all t G [0,1 ]\H*. 
Therefore the Borel set B = Uf^ 1 Bi has Lebesgue measure zero and for all i = 1,2,..., 


m = h(ty, vte[o,i]YB. 


Hence for every t G [0,1]\L>, 



x t ip(dx\t) = fi(t) 1 



x% d5 h(t), 


Vi = 1 , 2 ,.... 


where Sj ] o) is the Dirac measure at the point fift) G [0,1]. As measures on compact sets are 
moment determinate, one must have fj(dx\t) = for all t G [0,1 }\B. Therefore dnfx,t) = 

d f!(t)dt, i.e., the measure p is supported on the trajectory {(t,x(t)): t G [0,1]}, where xft) = 
fift) for almost all t G [0,1]. ■ 
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Remark 3.2. If the trajectory t t —> x(t) is a polynomial of degree say d, then the vector of 
Legendre-Fourier coefficients x E l 2 has at most d + 1 non-zero elements. Therefore for every 
j = 2,.. Xj E also has at most jd + 1 non-zero elements and the condition (3.5) can be 
checked easily. 


In Theorem 3.1(a) one assume that sup 

n 


£ * 0 ')-^ 


j=0 


< oo which is much weaker than, 

OO 


n 


e.g., assuming the uniform convergence 


£ j - x 
3=0 


—> 0 as n —> oo. The latter (which is 

OO 


also much stronger than the a.e. pointwise convergence) can be obtained if the function x(t) has 
some smoothness properties. For instance if x belongs to some Lipschitz class of order larger 
then or equal to 1/2, then uniform convergence takes place and one may even obtain rates of 
convergence; see, e.g., [12] and also [15] for a comparison (in terms of convergence) of Legendre 
and Chebyshev expansions. In fact, quoting the authors of [5], “... knowledge of the partial 
spectral sum of an L- 2 function in [—1,1] furnishes enough information such that an exponential 
convergent approximation can be constructed in any subinterval in which f is analytic”. 


Example 3.3. To illustrate Theorem 3.1 consider the following toy example with p on [0, l] 2 
and with marginal w.r.t. “t” being the uniform distribution on [0,1] and conditional f}{dx\t) = 
^exp(-t) f° r t £ [0) !]• That is, 1 i—» x(t) = exp(—t). 

Then the first 11 Legendre-Fourier coefficients x(j), j = 0,..., 10 of x read 


x = [0.63212055 -0.1795068 0.0230105 -0.0019370 0.0001217 

-0.0000061 0.0000002 -0.00000001 -0.00000004 -0.0000015 -0.0000625], 

Similarly the first 11 Legendre-Fourier coefficients of i i—>• x(t) 2 = exp(— 2t) read 


x 2 = [0.4323323 -0.2344075 0.0588678 -0.0097965 0.0012219 

-0.0001219 0.0000101 -0.0000007 0.00000004 -0.000000004 -0.00000004], 


With n 


5 the polynomial t e -> x 


( 2 ) 

5 


(t) := 




k =0 


reads 


10 

x 5 2) (*) = *5 2) ^) &3 (*) > teR, with 

k =0 

x = [0.4323336 -0.2344129 0.0588626 -0.0097976 0.0012219 

-0.0001218 0.0000098 -0.0000006 0.00000003 -0.000000001 0.0000000], 

and we can observe that x ^ — x 2 ~ O(l0 -5 ). 

In fact the curves i el a;) '(t) and t i —> exp(—2 1) are almost indistinguishable on the inter¬ 
val [0,1]. 


3.1 A more general case 

We have considered a measure p on [0, l] 2 whose marginal with respect to t E [0,1] is the 
Lebesgue measure. The conditions of Theorem 3.1 are naturally stated in terms of the (shifted) 
Legendre-Fourier coefficients associated with the functions t ^ f t (t) of L 2 QO, 1]) defined 
in (3.1). 

However, the same conclusions also hold if the marginal of p with respect to t E [0,1] is some 
measure dv = h(t)dt for some nonnegative function h E Li([0,1]) with all moments finite. The 
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only change is that now we have to consider the orthonormal polynomials 1 H y j = 0,1, 

with respect to v. Recall that all the 'Hj ’s can be computed from the moments 

7o O') = [ t j dfi(x , t) = [ t j du(t ) = [ t j h(t ) dt, j = 0 , 1 ,.... 

J Jo Jo 

Then proceeding as before, for every i = 1,2,..., 

[ V x l dfx(x,t) = f tJ ( [ x t, 4>(dx\t)]h(t) dt, j = 0 , 1 ,..., 

J Jo I 7 [o,i] J 

'-V-' 


and we now consider the vector of coefficients f hi = ( fhi(j )) defined by 

fhJj) = [ 'Hj{t)f i (t)h(t)dt, j = 0,1,... 

Jo 

the analogues for the measure du = h(t)dt and the function fi £ ([0, l],z/), of the (shifted) 
Legendre-Fourier coefficients fi(j) of fi in (3.1) for the Lebesgue measure on [0,1]. Then the 
conditions in Theorem 3.1 would be exactly the same as before, excepted that now, x = ( x(j )) 
with 


x(j) = [ Uj(t)x(t)h(t) dt, j = 0,1,.... 

Jo 

3.2 Discussion 

Theorem 3.1 may have some practical implications. For instance consider the weak formula¬ 
tion V of an optimal control problem P as an infinite-dimensional optimization problem on an 
appropriate space of (occupation) measures, as described, e.g., in [14]. In [8] the authors propose 
to solve a hierarchy of semidefinite relaxations (Vk), k = 1,2,..., of V. Each optimal solution 
of Vk provides with a finite sequence z fc = a g) such that when k —> oo, z k —> z* where z* is 
the infinite sequence of some measure d/r(t,x, u) on [0,1] x X x U, where X C M n , U C M m , 
are compact sets. 

Under some conditions both problems P and its relaxation V have same optimal value. If /a 
is supported on feasible trajectories {(t,x(t), u(t)) : t £ [0,1]} then these trajectories are optimal 
for the initial optimal control problem P. So it is highly desirable to check whether indeed fi is 
supported on trajectories from the only knowledge of its moments z* = (z* a p). 

By construction of the moment sequences z k one already knows that the marginal of fi with 
respect to the variable “t” is the Lebesgue measure on [0,1]. Therefore we are typically in the 
situation described in the present paper. Indeed to check whether /a is supported on trajectories 
{(t, (xi(t),..., x n (t), u\(t),..., u m (t)): t £ [0,1]}, one considers each coordinate Xi(t) or Uj(t) 
separately. For instance, for Xi(t) one considers the subset of moments 7 k(j) = i z j a 0 ) w ith 
j = 0,1,..., a = (0,..., 0, k, 0,..., 0) £ N n , k = 0,1,..., with k in position i. If (3.4) holds 
then indeed the marginal nt, Xi of fi on ( t,Xi ), with moments ( 7 k(j)) is supported on a trajectory 
{(t,Xi{t)): t £ [0,1]}. 

Of course, in (3.4) there are countably many conditions to check whereas in principle only 
finitely many moments of z* are available (and with some inaccuracy due to (a) solving numer¬ 
ically a truncation Vk of V, and (b) the convergence z k —> z* has not taken place yet). So an 
issue of future investigation is to provide necessary (or sufficient?) conditions based only on 
finitely many (approximate) moments of /a. 
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